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Abstract. - We consider the overdamped motion of Brownian particles, interacting via particle 
exclusion, in an external potential that varies with time and space. We show that periodic po- 
tentials that maintain specific position-dependent phase relations generate time-averaged directed 
current of particles. We obtain analytic results for a lattice version of the model using a recently 
developed perturbative approach. Many interesting features like particle-hole symmetry, current 
reversal with changing density, and system-size dependence of current are obtained. We propose 
possible experiments to test our predictions. 
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Q , Introduction. — Stochastic pumps refer to systems 
' Where overall directed motion of particles are obtained 
I ijnder the influence of external driving forces that are un- 

^ biased in the sense that they vanish either on spatial or 
temporal averaging. However, in order to generate di- 
rected motion, it is necessary that the external forces have 
some form of time-rcvcrsal-symmetry breaking inbuilt into 
.them. Natural examples of stochastic pump include ion- 
pumps associated with the cell-membrane (e.g., Na"*", K+- 
^^ATPase pump [1]), molecular motors (kinesin/dynein, 
I rnyosin) moving on polymeric tracks (microtubule, F- 

r*actin) [2], etc. Ref. [3] showed that the classical pump 
models for particle and heat transport [2,4-6] are similar 
^^to isothermal rachet models of molecular motors [2,7]. 

^ J Most studies of particle pumps focussed on systems of 
■ - non-interacting particles, though, there have been some 
numerical [8] and analytical work [9] on interacting Brow- 
nian motors. Recently a model of a classical pump, similar 
to those used in the study of quantum pumps [10, 11], has 
been proposed in Ref. [12,13]. They studied directed cur- 
rent (DC) in the presence of inter-particle interactions. 
The model studied was the symmetric exclusion process 
on a ring which closely mimics a system of diffusing parti- 
cles with short ranged repulsive interaction in one dimen- 
sion. It was shown that for time oscillatory hopping rates 
between two neighbouring points a DC current was estab- 
lished in the system. The hopping rates were taken to be 
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symmetric around every point, and this made it difficult 
to make a direct connection of this model to a system of 
particles in an oscillatory external potentials. In this letter 
we consider a stochastic exclusion dynamics which directly 
mimics hard core particles moving in a time-dependent ex- 
ternal potential. We show that the perturbative approach 
of Ref. [12, 13] can be used for this case as well, and thence 
obtain a number of interesting predictions for this system. 

Model. — The dynamics of small thermally diffusing 
interacting particles confined within a narrow tube is well- 
described by over-damped Langevin equations in one di- 
mension. Consider N particles interacting with each other 
through a short ranged repulsive interaction potential U 
and in an external potential V. Let Xr denote the po- 
sition of the r-th particle. The equations of motion are 
then given by: jdxr/dt = —dU /dxr — dV /dxr + , = 
1,2, ...,A , where £,r(t) is white Gaussian noise with 
= 0, (Cr(i)6(*')) = 2-ikBT5r,s5{t - t'), kB IS the 
Boltzmann constant and T the ambient temperature. We 
now discuss how this maps to a exclusion process. In the 
absence of an external potential, hard core colloidal par- 
ticles of diameter a, confined within a narrow channel of 
width < 2a, are well-described by the symmetric exclusion 
process (SEP) [14]. Discretizing space by the particle- 
size a, the SEP describes this system with a particle- 
hopping rate fo = D/a^ where the diffusion constant D 
obeys the Einstein relation D = ksT/^. The effect of 
an external potential V{x, t) is to make the hopping rates 
time- and position- dependent. Let P{x + + 5t\x,t) 



be the conditional probability density of a particle being 
at a; + a at time t + 6t, given that it is at x at time t. 
Then from the Langevin dynamics we see that the tran- 
sition rate is given by P{x + a,t + St \ x,t) a /St ^ 
fo[l - (Vix + a) - V{x))/2kBT], for weak potential differ- 
ences [(yix + a) - V{x)) < 2kBT]. 

Thus we consider a lattice model consisting of parti- 
cles on a ring of L sites with lattice spacing a. Each site 
I = 1, 2, . . . , i may contain ni = 0, 1 particle. A parti- 
cle at site I can hop to an empty nearest neighbour site 
I ± 1 with rate w/,i±i = /o[l — A(m;±i — u/)/2], where 
Am; = Vi/ksT. We consider driving at a constant fre- 
quency and set ui = ai sm{nt + (pi) = 2Re [rjie'^^*] where 
7]i = aie^'^'-/2i. For time-independent potential Vi, i.e., 

= 0, the transition rates satisfy detailed balance and 
the system reaches thermal equilibrium with zero current. 
As we show now, in the presence of time-dependent V;, 
the system is driven out of equilibrium and can show di- 
rected flow of particles resulting from AC driving forces. 
We analyze the model utilizing a pcrturbative expansion 
in the strength A of time-periodic potential, and also via 
Monte-Carlo simulations. 

Calculation of the current. The current from 
site I — 1 to I, Ji-i^i is defined by the continuity equation 
obtained from the time evolution of local density pi = 
(ni), where (...) denotes an average over the noise. This 
gives [15] dpi/dt = Ji^ij - Ji^i+i where 

J(_l.; = (Wi-l,i Pl-l - Wi,i-i pi) - {Wi-i,i - Wi,i-i) Ci-i^i 

(1) 

gives the noise averaged current from site I — 1 io site 

1 and Ci^m = {n-irim) denotes the two-point correlations. 
Ci,m appears in the expression of Ji-ij due to the particle 
exclusion. In the time-periodic steady state the current 
on each bond, averaged over the time period r = 27r/f2, 
is the same and hence the net DC current is given by 
J = {1/Lt)J2i'=i Jq dtJi-i^i{t). Using the form of the 
hopping rates we then get: 

J^-^y^f dtiui~ui^i)ipi^i+pi-2Q-ij). (2) 
1=1 -l^ 

To evaluate J we need to evaluate pis and C;_i^;s and 
this we do using perturbation theory. For weak driving 
potential we make a perturbative expansion of pi, Ci^m in 
a series in the dimensionless parameter A 
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k=1.2, 



The A = terms in the above expansions correspond to 
the absence of any external potential and this is then 
just the symmetric exclusion process with uniform hop- 
ping rate /q. For that case the steady state is the equilib- 
rium state obeying detailed balance and is characterized 



by the position-independent densities and correlations. 
These are all known exactly, for example p = N/ L = p, 
cl% = p{N-l)/iL-l),etc. [15]. 

As noted in Ref. [12] the time evolution of the first or- 
der terms p\^\t) and cl^i(t) are given by the following 
equations: 

see eq. (4), (5) 

where Aigi^^n = 9i+i.m + gi-i,m-2gi^m defines the discrete 
Laplacian, and qo = p-C[^^, ko = C'[°^-C[°^ ^ with equi- 
librium three-point correlation C;^*^ „ = cl^^{N — 2)/{L — 
2). These linear homogeneous equations can be solved ex- 
actly to obtain the long-time oscillatory solution [12]. We 
find: 

pp'(i) = 2Re[^ze'"*] (7) 
where the vector A = {Ai, A2, . . . , Al} is given by 



A = 



qofo 



in - /oA 



iilqQ 



in - foA 



V 



(8) 



with T] = {rji, 772, . . . , and A the discrete Laplacian op- 
erator. The eigenfunctions of A are given by (l/-\/L)e^"'' 
with corresponding eigenvalues eg = —2(1 — cosq). Ex- 
panding A using its cigcnbasis we get Ai — --qorji + 

Em^hmVm whcrc Rl^„^ = (^OgoA) Eg e*'^'""V(*^^ - 
fotq) with q = 27rjyL, j = 1,2, . . . , L. The equation of 
two-point correlations can also be solved and it is easy to 
verify the following steady state solutions (for all l,m): 



clUit) = %['\t) + p^^Ht)]. 

Qo 



(9) 



Then using Eq. 2 one can find the general expression for 
the directed current. 
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where fjg = (l/y/L) E 

Note that the effect of interaction is entirely contained 
in the prefactor (go — 2fco), which in the large L limit equals 
to p(l — p) (1 — 2p) . Several interesting features follow. The 
current vanishes at half-filling and its sign reverses with in- 
creasing particle density. The dynamics has particle-hole 
symmetry, and this is explicit in the density dependence. 
At low densities we obtain J cx p and this corresponds to 
the case where interactions can be neglected. This differs 
from the result in Ref. [12] where the current vanished 
in the absence of interactions. We now investigate the 
solution in Eq. 10 for different choices of the oscillating 
potential. We consider the following two cases: 



dp, 
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/o(A; + A^)CIII + fokoiAiui + AmUm) for / 7^ m ± 1, 



dt 



dC, 
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da 



dt 
(1) 



,1+1 



/o(C/-i,/+i + C'u+2 - '^^Li'+i) + foko{ui-i + Ui+2 -ui- Ui+i) 
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dt 
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Fig. 1: (Color online) Directed current J2 as a function of mean 
density p. The points denote Monte-Carlo results and the line 
is a plot of the function in Eq.ll. The parameters used are: 
system size L = 16, hopping rate fo = 0.34, potential strength 
A = 0.5, time period of oscillation r = 10 (il = 0.2tt), and 
phase difference (j> = — 7i"/2. The data (points) were collected 
over lOOr after equilibration over another lOOr. All the data 
were averaged over 10* initial conditions. 



Case (i) : Localized pump. Here we consider 
the specific case of a localized pump with time-varying 
potentials acting only on two consecutive sites such that 
ai = tti = 1 with all other ai = 0, and (fii — (f), (pL — 0- 
For this case Eq. 10 leads to the DC current 



„ , .fl sin 6 ^ 



(11) 



where_z_ = ?//2 - ^/[y/2f - 1 with y^2 + in/fo. Note 
that J2 has a sinusoidal dependence on the phase differ- 
ence of external driving 0, leading to maximal driven cur- 
rent a.t (f> = ±7r/2. J2 decays as 1/L with the system size 
L. 

In Fig. 1 we plot density dependence of J2 obtained from 
Monte-Carlo simulations. The data show good agreement 
with Eq. 11. At p = 1 the current vanishes due to com- 
plete jamming of particles. Note the zero in the current at 
the point of half- filling p = 1/2, and the reversal in the di- 
rection of current as the density crosses this point (Fig. 1). 
The maximal currents are obtained at p = (1 ± l/V3)/2. 
We re-emphasize that the density dependence is general, 
and independent of whether the pumping is localized or 
spatially distributed. 



Fig. 2: (Color online) System-size independent currents J2L 
and Jl as a function of frequency Q using Eqs. 11 and 12 with 
L — 10* and p — {I ~ l/\^)/2. All other parameters are the 
same as in Fig. 1. The maximum in Jl is at = 2/o — 0.68. 



Case (ii) : Traveling wave. — A spatially 
distributed external potential of the form V{x,t) = 
Vq sm{^lt + (j)x) where (j) is constant generates a directed 
traveling wave of external force. We use ai ~ 1 and 

= with a constant 0, i.e., rji — exp(i0Z)/2i for all I. 
Hence we have fjq = {l/2i)^/LSq^^, and using this expres- 
sion in Eq. 10 we find the time-averaged current 



Jl 
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= -A2(go-2fco)/o' 



il sin 0(1 — cos 
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(12) 



Since in this case external force acts on all the sites of the 
system Jl is independent of system size L. This behaviour 
is in contrast to the decay of current in localized pump as 
J2 (X 1/L. While Eqs. 11 and 12 contain the full func- 
tional dependence of the time-averaged currents for the 
two cases, it is instructive to illustrate their frequency- 
dependence at a density and phase at which the currents 
are at their maximum. Thus, using (p — —tt/2 for a sys- 
tem of density p = (1 — 1/^/3)/ 2 and a large system size 
L ~ 10"', the dependence of J2i, Jl on driving frequency 
is shown in Fig. 2. At small frequencies, current grows 
linearly with fi. Eq. 12 shows that Jl reaches a maximum 
value J™'''' = A2/o/24V3 at a driving frequency = 2/o. 
This denotes a resonance between the intrinsic relaxation 
rate /o and the external driving frequency £7. At high 



frequencies ft > 2/o, the relaxation of individual parti- 
cles lags behind the fast variation of external force, and 
eventually the current decays as 1/0. 

Outlook. — Our model of stochastic exclusion dynam- 
ics showed emergence of DC current in response to peri- 
odic forcing that vanishes on temporal averaging. Definite 
phase-difference between the oscillatory forcing at neigh- 
bouring sites breaks the time-reversal symmetry to yield 
directed current. The mechanism works even in the limit 
of very low densities where particles rarely interact. How- 
ever, with increasing density we found a non-nionotonic 
variation of the DC current, entirely due to the exclu- 
sion interaction between particles. Most interestingly, the 
system showed a current reversal around half filling: At 
density > 1/2, the driven DC current flows in a direc- 
tion opposite to that in the non-interacting limit. This 
behavior is associated with the particle-hole symmetry - 
the dynamics of holes at high densities is equivalent to the 
dynamics of particles at low densities. 

Our predictions may be tested experimentally by us- 
ing oscillatory force on colloids moving in confined ge- 
ometries. There have been several experimental obser- 
vations in colloidal systems of single-file-difFusion which 
is one of the signatures of diffusing interacting particles 
moving in one dimension [14,16-19]. It is possible to ap- 
ply localized AC drive on such confined-colloids by, e.g., 
exploiting AC electro-kinetic effects, or localized oscil- 
latory pressure. Alternatively in the setup of Ref. [18] 
one could impose oscillatory trapping potential by using 
laser traps of tunable power. Consider the motion of 
colloidal particles of diameter a « 2/^m moving in wa- 
ter inside a narrow tube of diameter < 2a. At room 
temperature ksT = 4.2 x 10"^"'^ N-m using water viscos- 
ity ly « 10~^ Ns/m^ one finds the Stokes-Einstein self- 
diffusion constant Dq = ksT /'inva « 0.2/im^/s. This 
gives hopping rate /o = Do/a^ = 0.05 s^^. Then for the 
case of extended driving corresponding to the result in 
Eq. 12, the maximum current is obtained at a density 
of p = (1 — l/\/3)/2 w 0.2 and a driving frequency of 
il = 2/o w 0.1 Hz. Since the value of the maximal current 
Jl ~ A'^/o/40, this means that we can get particle flow 
velocities ^ J^a ~ 0.0025 fim/a. This compares with the 
drift velocity that would be attained by a colloidal particle 
of the same radius and carrying a charge of 10 electrons, 
when placed in an electric field ^ 30V/m. 
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